Boundary effects on energy dissipation in a cellular automaton model 
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In this paper, we numerically study energy dissipation caused by traffic in the Nagel-Schreckenberg 
(NaSch) model with open boundary conditions (OBC). Numerical results show that there is a non- 
vanishing energy dissipation rate Ed, and no true free-flow phase exists in the deterministic and 
nondeterministic NaSch models with OBC. In the deterministic case, there is a critical value of 
the extinction rate f3 c d below which Ed increases with increasing /3, but above which Ed abruptly 
decreases in the case of the speed limit w max ^ 3. However, when w max ^ 2, no discontiguous 
change in Ed occurs. In the nondeterministic case, the dissipated energy has two different contribu- 
tions: one coming from the randomization, and one from the interactions, which is the only reason 
for dissipating energy in the deterministic case. The relative contributions of the two dissipation 
mechanisms are presented in the stochastic NaSch model with OBC. Energy dissipation rate Ed is 
directly related to traffic phase. Theoretical analyses give an agreement with numerical results in 
three phases (low-density, high-density and maximum current phase) for the case v max = 1 • 

PACS numbers: 05.65.+b, 45.70.Vn, 05.60.-k, 89.40.Bb 



I. INTRODUCTION 



In the last decades, traffic problems have attracted 
much attention of a community of physicists because of 
the observed nonequilibrium phase transitions and vari- 
ous nonlinear dynamical phenomena. A number of traffic 
models have been proposed to investigate the dynami- 
cal behavior of the traffic flow, including fluid dynami- 
cal models, gas-kinetic models, car-following models and 
cellular automata (CA) models 0, H, & Q. These dy- 
namical approaches represented complex physical phe- 
nomena of traffic flow among which are hysteresis, syn- 
chronization, wide moving jams, and phase transitions, 
etc. Among these models, the cellular automata ap- 
proaches can be used very efficiently for computers to 
perform simulation % USB S I, M El Hi El El • 
The Nagel-Schreckenberg (NaSch) model is a basic CA 
models describing one-lane traffic ffowQ. Based on 
the NaSch model, many CA models have succeeded 
in modeling a wide variety of properties of vehicular 
traffic 0, 1, 0, 1, i M EL El El, EI- 

On the other hand, the problems of traffic jams, en- 
vironmental pollution and energy dissipation caused by 
traffic have become more and more significant in mod- 
ern society. Financial damage from traffic due to en- 
ergy dissipation and environmental pollution is huge ev- 
ery year. In accordance with previously reported re- 
sults in Refs.[15], more than 20% fuel consumption and 
air pollution is caused by impeded and "go and stop" 
traffic. Most recently, the problem of energy dissi- 
pation in traffic system has been investigated in the 
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framework of car following model, city traffic model and 
NaSch model with periodic boundary conditions (PBC), 
respectively [H ES El, El H3]- And analytical expres- 
sions for energy dissipation have also been provided in 
the nondeterministic NaSch model with PBC in the case 
of v max = 1 and free-flow statc[20]. However, the ef- 
fects of boundary condition on energy dissipation have 
not been discussed yet. 

The most significant difference between systems with 
open and periodic boundary conditions is the vehicle 
density p. In a periodic system, which has no maxi- 
mum current phase, vehicle density is considered as an 
adjustable parameter. In systems with open boundary 
conditions (OBC), however, there are two adjustable pa- 
rameters, namely the injection rate a and the extinc- 
tion rate (5, and the vehicle density p is only a derived 
parameter. Compared with periodic systems, it implies 
that open systems show a different behaviour of quan- 
tities such as the global density, the current, the den- 
sity profile and even the microscopic structure of traffic 
phase[H [H [H [H [H [H, [H M, [H, [H. Therefore, 
energy dissipation in the CA model with OBC, which is 
relevant to many realistic situation in traffic, should be 
further investigated. 

In this paper, we investigate the energy dissipation rate 
within the framework of the deterministic and nondeter- 
ministic NaSch model with OBC. The behaviours of the 
energy dissipation rate in different traffic phase are dis- 
tinct. Theoretical analyses are presented in low-density, 
high-density and maximum current phase in the case of 
Vmax = f ■ The influences of the speed limit i> max on the 
energy dissipation are also investigated. Energy dissipa- 
tion caused by braking is related not only to the velocity 
of vehicles, but also to the headway distribution. Thus, 
the behaviour of the energy dissipation rate is more com- 
plex than simpler quantities, and should be further inves- 
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tigated. The results of this article may lead to a profound 
understanding of some features of traffic system or may 
provide schemes for reducing energy dissipation of the 
existing traffic network. 

The paper is organized as follows. Section II is de- 
voted to the description of the model and the definition 
of energy dissipation rate. In section III, the numerical 
studies are given, and the influences of the injection and 
extinction rate on energy dissipation rate are considered. 
And theoretical analyses are presented in the special case 
"Umax = 1. Finally, the conclusions are given in section 
IV. 



II. MODEL AND ENERGY DISSIPATION 

Our investigations are based on a one dimensional cel- 
lular automaton model introduced by Nagel and Schreck- 
enberg. The model is defined on a single lane road con- 
sisting of L cells of equal size numbered by i = 1, 2, • • • , 
L and the time is discrete. Each site can be either empty 
or occupied by a car with the speed v = 0, 1, 2, • • ■ , u max , 
where w m ax is the speed limit. Let x(i,t) and v(i,t) de- 
note the position and the velocity of the ith car at time 
t, respectively. The number of empty cells in front of the 
ith vehicle is denoted by d(i, t) — x(i + 1, t) — x(i, t) — 1. 
The following four steps for all cars update in parallel 
with periodic boundary. 

(1) Acceleration: 

u(i,t + l/3) -> vaia[v(i,t) + l,w max ]; 

(2) Slowing down: 

v(i,t + 2/3) -> xaux[v(i,t+ 1/3), d(«,t)]; 

(3) Stochastic braking: 

v(i, t+1) — > max[v(i, t+2/3) — 1, 0] with the probability 

p; 

(4) Movement: x(i, t + 1) — > x(i, t) + v(i, t + 1). 
Open systems are characterized by the injection rate 

a and the extinction rate /3, which means by the proba- 
bility a and (5 that a vehicle moves into and out of the 
system. In this paper, open boundary conditions are de- 
fined according to [28|, [29|. At site i = which means out 
of the system, a vehicle with speed v = v max is created 
with probability a. The vehicle immediately moves for- 
ward in accordance with the NaSch rule. If the site i = 1 
is occupied by a car, the injected vehicle at site i = is 
deleted. At i — L + 1 a "block" occurs with probability 
1 — j3 and causes a slowing down of the vehicles at the 
end of the system. Otherwise, a vehicle may leave freely 
from the end of the system. 

It should be mentioned that for i> max = 1 the model 
above in the nondeterministic case is different from par- 
allel updated asymmetric exclusion pro cess (ASEP) with 
open boundary conditions pU. I22I l23j|. In the ASEP 
model, if the site L is occupied then the particle on that 
site exits with probability /3, irrespective of their veloc- 
ity. In the model above, however, the extinction of the 
vehicle at site L depends not only on the probability /3, 
but also on the braking probability p. Even if the vehi- 



cle with v = 1 is at site L and there is no "block" at 
site L + 1, it may fail to exit because of being randomly 
delayed. This difference may influence the current and 
phase transition, which will be analysed in the following 
section. In the deterministic case, however, the model 
above with v max — 1 is identical with parallel updated 
ASEP with OBC, for there is no stochastic delay. 

The kinetic energy of the vehicle with the velocity v is 
mv 2 /2, where m is the mass of the vehicle. When braking 
the energy is lost. Let Ed denotes energy dissipation rate 
per time step per vehicle. For simple, we neglect rolling 
and air drag dissipation and other dissipation such as 
the energy needed to keep the motor running while the 
vehicle is standing in our analysis, i.e., we only consider 
the energy lost caused by speed-down. The dissipated 
energy of ith vehicle from time t — 1 to t is defined bv[2p| 

_ ff [v 2 {i,t- 1) -v 2 {i,t)] for v{i,t) < v(i,t - 1) 
6(Z ' ' ~ \ for v(i,t) ^ v(i,t- 1). 

Thus, the energy dissipation rate 

t +T N 

E *=TN £ 5> (M) ' (2) 
t=t +i t=i 

where N is the number of vehicles in the system and 
<o is the relaxation time, taken as to = 10 5 . In this 
model, the particles are " self-driven" and the kinetic en- 
ergy increases in the acceleration step. In the stationary 
state, the value of the increased energy while accelerat- 
ing is equivalent to that of the dissipated energy caused 
by speed-down, and the kinetic energy is constant in the 
system. In the simulation, the system size L = 1000 is se- 
lected, and the results are obtained by averaging over 20 
initial configurations and 10 4 time steps after discarding 
10 5 initial transient states. 



III. NUMERICAL RESULTS 

A. Effects of the boundary conditions on energy 
dissipation in the deterministic case 

First, we investigate the influences of the boundary 
conditions on energy dissipation in the deterministic 
NaSch model with the maximum velocity w max = 5. In 
the deterministic case, the stochastic braking is not con- 
sidered, i.e., p = 0. Figure 1 shows the energy dissipation 
rate Ed as a function of the extinction rate (3 with dif- 
ferent values of the injection rate a. As shown in Fig. 
1, there is a critical value of the extinction rate (3 c d be- 
low which Ed increases with the increase of the rate /3, 
but above which Ed abruptly decreases. The position of 
(3 c d shifts towards a high value of [3 with increasing the 
injection rate a. The astonishing result is that there is 
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FIG. 1: Energy dissipation rate E d (scaled by m ) as a func- 
tion of the extinction rate in the deterministic NaSch model 
with w max = 5 for various values of the injection rate a. 



FIG. 2: Energy dissipation rate E d (scaled by m ) as a func- 
tion of the extinction rate in the deterministic NaSch model 
in the case of a = 1.0 for various values of the speed limit 

Wmax - 



a nonvanishing energy dissipation, even though in low- 
density phase, which means that there is no "true" free- 
flow phase. 

When < cd , the state of the system is high-density 
phase. As the extinction rate increases, the kinetic 
energy possessed by vehicles increases because of the in- 
crease of the mean vehicle velocity, and the dissipated 
energy while braking increases. When > c d, the state 
of the system is low-density phase in which the distance- 
headways is larger and the interaction between vehicles is 
weaker than that in the high-density phase. The interac- 
tion is the only reason for dissipating energy in the deter- 
ministic case. Consequently, Ed decreases abruptly when 
the transition from high-density to low-density phase oc- 
curs. Because of boundary effects there are vehicular in- 
teractions even at low-density phase. Therefor, there is a 
nonvanishing energy dissipation rate Ed, and no "true" 
free-flow phase exists in the deterministic NaSch model 
with open boundary conditions. 

Energy dissipation rate Ed in the case of v max = 3 
and 4 show similar behaviour to that for v max = 5, and 
the position of c d shifts towards a low value of with 
increasing the speed limit u max , as shown in Fig. 2. But 
in the case of w max = 2 and 1, there is no critical value 
of the extinction rate c d, i.e., no discontinuous change 
in Ed occurs. According to previously reported results in 
Refs.[19], in the case of v max < 3, the state of a system 
with injection rate a = 1 is the jamming phase, while 
in the case of i> max ^ 3, the jamming state exists in 
the region of low value of 0, and low-density lies in the 
region of high value of 0. Thus, there is no discontinuous 
change in Ed in the case of w max = 2 and 1 for no phase 
transition occurs. Near = 0, energy dissipation rate Ed 
is independent of the speed limit u max , and Ed increases 
linearly with the increase of 0. In the case of a = 1 
and — > 0, the maximum velocity which vehicles can 
move is 1, so the speed limit has no influences on energy 
dissipation. 

In the special case of w max = 1, energy dissipation rate 
Ed is proportional to the mean density of " go and stop" 




FIG. 3: Energy dissipation rate Ed (scaled by m ) as a func- 
tion of the extinction rate in the deterministic NaSch model 
with w m ax = 1 for various values of the injection rate a. Sym- 
bol data are obtained from computer simulations, and solid 
line corresponds to analytic results of the formula (6). 

vehicles per time step. The mean density of "go and 
stop" vehicles is defined in the following way: 



1 1 



t +T N 



(3) 



t=t + l i=l 



where n(i,t) = for stopped cars and n(i,t) = 1 for 
moving cars at time t, and to is the relaxation time as 
mentioned in the section II. And energy dissipation rate 
Ed in the case of w max — 1 can be written as 



E d = -jmp gs . 



(4) 



For open boundary conditions, the mean density of "go 
and stop" vehicles per time step reads 



P gs = n - nl = 0(1 - 0), 



(5) 



where n = N /N is the fraction of the stopped vehicles, 
and A^o is the number of stopped vehicles on the road. 
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As a consequence, energy dissipation rate Ed in the 
case of » m „ — 1 can be obtained as 



(6) 



In formula (6), energy dissipation rate Ed is directly re- 
lated to the probability for a vehicle moving out of the 
system or the probability for a vehicle occupying the last 
site of the system. As shown in figure 3, theoretical anal- 
ysis is in good agreement with numerical results. 

For t> max = 1 and a < 1 , a discontinuous change in Ed 
also occurs at a critical point (3 c d, and the position of (3 c d 
shifts towards a high value of (3 with increasing a, which 
is similar to that in the case of u max = 5. 



B. Effects of the boundary conditions on energy 
dissipation in the nondeterministic case 

Next, we investigate the rate of energy dissipation Ed 
when the stochastic braking behaviours of drivers are 
considered, i.e., p ^ 0. In the nondeterministic case, the 
dissipated energy has two different contributions: one 
coming from the stochastic noise and one from the inter- 
actions between vehicles (this in fact is the only reason 
for dissipating energy in the deterministic case) . Let Edi 
and Ed r denote the rate of energy dissipation caused by 
the interactions and randomizations, respectively. And 
energy dissipation rate Ed = Edi + Ed r - Figure 4 and 5 
show the relation of the energy dissipation rate Edi and 
Edr to the extinction rate f3 with different values of the 
injection rate a, respectively, in the case of u max = 5 and 
p= 0.5. As shown in figure 4, there is a critical value of 
the extinction rate [3 cr below which Edi increases with 
increasing /3, but above which Edi abruptly decreases, 
except for the case of a = 1. Different from Edi, energy 
dissipation rate Ed r sharp increases at the critical point 
(3 cr above which Ed r shows the approximate plateau and 
is independent of (3, except for a = 1, as shown in figure 
5. And the position of (3 cr shifts towards a high value of 
[3 with the increase of the injection rate a. For a = 1, 
there is no discontinuous change in Edi and Ed r - Com- 
pared figure 4 with 5, it turns out that energy dissipation 
is mainly caused by the randomization in the low-density 
phase, and by the interactions in the high-density phase. 
When a > 0.35, the values of the rate of energy dissipa- 
tion Edi and Ed r for various value of a collapse into a 
single curve (not shown). 

In fact, in the case of p= 0.5 and v max = 5, there are 
the transitions from the high-density to low-density phase 
for a ^ 0.35 and from high-density to maximum current 
phase for a > 0.35[29j . In the high-density phase, with 
the increase of the extinction rate /?, the mean velocity 
increases and the dissipated energy increases. In the low- 
density and maximum current phase, with increasing the 
rate /?, the distance headways increases and the interac- 
tions lowers; thus energy dissipation rate Edi decreases. 
Energy dissipation rate Ed r reaches a constant value in 
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FIG. 4: Energy dissipation rate Edi (scaled by m ) as a func- 
tion of the extinction rate /3 in the non-deterministic NaSch 
model with u ma x = 5 and p = 0.5 for various values of the 
injection rate a. 




FIG. 5: Energy dissipation rate Edr (scaled by m ) as a func- 
tion of the extinction rate /3 in the non-deterministic NaSch 
model with w max = 5 and p — 0.5 for various values of the 
injection rate a. 



the low-density phase, but continuatively increases with 
the increase of fj in the maximum current phase. 

It should be noted that energy dissipation rate Ed for 
a = 1 is minimum, i.e., energy dissipation in the maxi- 
mum current phase is lower than that in the low-density 
phase (not shown). Traffic flow, however, in the maxi- 
mum current phase is maximal. 

The relationship of Ed to the extinction rate (3 for dif- 
ferent values of D max in the case of a = 1 is shown in figure 
6. Similar to the deterministic case, when the value of 
(3 is very small, the rate of energy dissipation shows a 
scaling relation and is independent of the speed limit, as 
shown in figure 6. In the region of high values of /?, the 
scaling relations of energy dissipation Ed to the speed 
limit cannot be observed; and Ed increases with the in- 
crease of v max . When w max > 3, energy dissipation rate 
Ed increases with increasing the rate (3 in the region of 



high values of (3. However, when 



^ 3, the value of 



Ed tends to be invariable. Though the states of the sys- 
tem for different speed limit are maximum current phase 
in the region of high values of (3, with the increase of (3, 



FIG. 6: Energy dissipation rate E d (scaled by m ) as a func- 
tion of the extinction rate j3 in the non-deterministic NaSch 
model with a — 1.0 and p — 0.5 for various values of the 
speed limit « max . 

the mean velocity increases for w max > 3 but does not 
vary in the case of v max ^ 3. Consequently, there is a 
plateau for the case v ma x ^ 3 in the region of high values 
of /?, which is different from that in the case of w max > 3. 

In the system with w max = 1, some vehicles can be 
stopped due to the stochastic braking, therefore the rate 
of energy dissipation Ed is proportional to the mean " go 
and stop" density p gs , which demonstrates that the prob- 
ability for "go and stop" vehicle to appear per time step 
in the system. 

In the low-density phase A (a < (5, a < a c = 1 — ^/p), 
the fraction of the stopped vehicles reads 



where q = 1 — p. And the mean " go and stop" density 
p gs can be obtain as 

2 Q-ot A q-a\ 

Pgs = n Q - n Q = 1 - . (8) 

1 — a \ 1 — a J 

Substituting formula (8) into (4), we can obtain energy 
dissipation rate Ed in the low-density phase A 

A m{q-a){l-q) 
Ed= 2(1 - a )2 ■ (9) 

Figure 7 shows the relation between the rate of en- 
ergy dissipation Ed and the injection rate a with various 
values of p, in the case of f3 = 1. As shown in Fig. 7 
in the region of low values of a, formula (9) gives good 
agreement with the simulation data. 

From -J^- = 0, we obtain the critical value a mc = 
1 — 2p in which the curve described by formula (9) 
reaches the maximum. Compared with the critical rate 
a c = 1 — y/p, the critical value of the stochastic braking 
probability Pc — j can be obtained. When p < p%, en- 
ergy dissipation rate E^ increases with increasing the in- 
jection rate a. When p ^ 0.5, however, with the increase 



FIG. 7: Energy dissipation rate E d (scaled by m ) as a func- 
tion of the injection rate a in the non-deterministic NaSch 
model in the case of i; max = 1 and (3 — 1.0 for various values 
of the stochastic braking probability p. Symbol data are ob- 
tained from computer simulations, and solid line corresponds 
to analytic results of the formula (9) and (15). 
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FIG. 8: Energy dissipation rate Ed (scaled by m ) as a func- 
tion of the extinction rate (3 in the non-deterministic NaSch 
model in the case of i) m « = 1 and a — 1.0 for various values 
of the stochastic braking probability p. Symbol data are ob- 
tained from computer simulations, and solid line corresponds 
to analytic results of the formula (12) and (15). 

of the rate a, energy dissipation rate E^ decreases. In 
the interval 0.5 > p ^ p^ , with increasing a, the rate of 
energy dissipation E^ increases first and decreases after 
a maximum value is reached. 

In the high-density phase B (/3 < a, f3 < [3 C — 

the model of this paper is different from the ASEP with 
parallel update and the mean velocity is determined not 
only by the extinction rate (3, but also the stochastic 
braking probability p. The fraction of the stopped vehi- 
cles reads 

n = 1 - qP. (10) 

And the mean "go and stop" density p gs can be written 
as 

Pgs =n -n 2 = q(3(l - q(3). (11) 
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Substituting formula (11) into (4), we can obtain the rate 
of energy dissipation Ed in the high-density phase B 

Ei = ^(q0-q 2 n (12) 

Figure 8 shows the relation the energy dissipation rate 
Ed as a function of the extinction rate /3 with different 
values of the stochastic braking probability p, in the case 
of a = 1. As shown in Fig. 8, the agreement can be 
obtained in the case of low value of the extinction rate (3. 

From -J^- = 0, we can obtain the critical value (3 mc = 
2(i^ p ) m which the curve corresponding to formula (12) 
reaches the maximum. Compared with the critical rate 
f3 c = the critical value of the stochastic braking 

probability p^ — \ can be obtained. When p ^ pjf, en- 
ergy dissipation rate E^ increases with the increase of 
the extinction rate j3. When p < p^, however, with in- 
creasing the rate (3, energy dissipation rate Ef increases 
first and decreases after a maximum value is reached. 

In the maximum current phase C (a > a c ,(3 > (3 C ), 
the fraction of the stopped vehicles reads 

n = y/p. (13) 

And the mean "go and stop" density p gs can be obtain 
as 

P gs = n Q - nl = y/p{l - y/p). (14) 

Substituting formula (14) into (4), we can obtain energy 
dissipation rate Ed in the maximum current phase C 

E c d =^(^-p). (15) 

Equation (15) demonstrates that energy dissipation in 
the maximum current phase is independent of the rate a 
and /3, and is only determined by the stochastic braking 
probability. Figure 7 and 8 give a comparison between 
number results and Eq. 15. As shown in the right region 
of figure 7 and 8, formula (15) gives good agreement with 
the simulation data. 



IV. SUMMARY 



In this paper, we investigate the rate of energy dissi- 
pation caused by braking in the NaSch model with open 
boundary conditions. Different from periodic systems, 
open systems in which the vehicle density is only a de- 
rived parameter arc controlled by the injection and ex- 
tinction rate. In fact, real traffic systems are usually 
open, hence it is highly desirable to investigate energy 
dissipation in traffic systems both numerically and theo- 
retically. 

Numerical results show that in the deterministic case 
there is a critical value of the extinction rate (3 c d above 
which Ed decreases abruptly for w max > 3, however, no 
discontiguous change in Ed occurs when v max < 3. The 
rate (3 c d is related not only to the injection rate a, but 
also to the maximum velocity of vehicles. In the nondc- 
terministic case, there is also a critical value of the ex- 
tinction rate (3 cr below which Edi and Ed r increase with 
increasing /?, above which Edi abruptly decreases but Ed r 
sharp increases and shows the approximate plateau with 
further increase of /3, when the transition from the high- 
density to low-density phase occurs. However, when the 
transition from the high-density to the maximum cur- 
rent phase occurs, the values of energy dissipation rate 
Edi and Ed r for various value of a collapse into a single 
curve, and no discontinuous change occurs. Moreover, 
there is a nonvanishing energy dissipation rate, and no 
" true" free-flow phase exists in the deterministic and non- 
deterministic NaSch models with open boundary condi- 
tions. 

Energy dissipation rate Ed is directly related to traf- 
fic phase. Energy dissipation in maximum current phase 
is smaller than that in the low-density phase. A phe- 
nomenological mean-field theory is presented to describe 
the energy dissipation rate Ed in three phases (low- 
density, high-density and maximum current phase) in the 
case of w max = 1. Theoretical analyses give an excellent 
agreement with numerical results. But in the case of 
u max > 1, explicit expressions about the energy dissipa- 
tion rate Ed do not be obtained because of effects of long 
length of time space correlations, and deserve further in- 
vestigate. 
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